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appear  in  J.  Opt.  Soc.  Am.  H  (1986).  Presented  at  the  Annual 
Meeting  of  the  Optical  Society  of  America  (Seattle,  WA,  October 
1986). 

2.  P.  L.  Marston  and  S.  C.  Billette,  "Scattering  of  light  by  a  coated 
bubble  in  water  near  the  critical  and  Brewster  scattering  angles" 
(abstract  only),  J.  Acoust.  Soc.  Am.  Suppl.  &Q,  59  (1986).  To  be 
presented  at  the  1 12th  Meeting  of  the  Acoustical  Society  of  America 
(Anaheim,  CA,  Dec.  1986). 

E.  Other  Written  Communications 

1 .  Informal  progress  letter  by  P.  L.  Marston  to  M.  Y.  Su  (dated  April  25, 
1986). 

2.  Letter  by  P.  L.  Marston  to  M.  Y.  Su  (dated  June  17, 1986)  on  the 
generation  of  microbubbles  for  calibration  purposes. 

II.  INTRODUCTION  AND  REVIEW 

Previous  research  (supported  by  the  Office  of  Naval  Research,  Physics  Division) 
considered  the  understanding  of  the  optical  scattering  properties  of  microbubbles  in 
water.1*7  The  present  research  involves  the  application  and  extension  of  that 
understanding,  to  recommendations  concerning  the  design  and  operation  of  an  optical 
bubble-size  spectrometer.  The  intended  instrument  would  measure  the  sizes  and  number 
densities  of  bubbles  in  ocean  waters.  The  data  are  to  be  recorded  real-time  and  it  would  be 
desirable  that  the  instrument  facilitate  remote  operation.  Hence  photography  of 
microbubble  populations  was  judged  to  be  unsuitable  whereas  the  electronic  recording  of 
optical  scattering  data  from  bubbles  should,  in  principle,  be  invertible  to  give  the  desired 
size  measurement.  The  resources  of  the  present  contract  ($10,000)  were  directed  entirely 


towards  theoretical  problems  of  interest  for  the  design  and  operation  of  such  a 
spectrometer. 

So  as  to  better  illustrate  the  contributions  of  the  present  research,  it  is  appropriate  to 
note  that  the  previous  research1'9  used  the  following  three-pronged  approach  to 
understanding  scattering  from  bubbles:  (i)  quantitative  and  qualitative  observations  and 
discovery  of  new  scattering  phenomena;  (ii)  numerical  evaluation  of  partial-wave  series 
(such  as  the  Mie  series)  which  are  for  special  bubble  shapes  (e.g.,  uncoated  spheres)  and 
comparison  of  these  results  with  data;  (iii)  construction  and  evaluation  of  simple  physical 
models  which  are  both  quantitatively  useful  and  promote  understanding.  The  general 
motivation  for  such  research  is  reviewed  briefly  in  Ref.  3, 4,  and  6.  Recent  experiments, 
carried  out  here  at  Washington  State  University,  have  emphasized  the  backscattering  of 
light  from  bubbles  in  water.6*7  Backscattering  may  be  particularly  useful  if  the  optical 
source  and  detector  can  not  be  placed  close  to  the  water  sample  which  contains  the  bubbles. 

For  the  instrument  design  under  consideration  by  NORDA,  however,  the  optical 
detectors  and  source  may  lie  within  a  few  feet  (or  closer)  of  the  ocean  water  sample.  It 
appears  that  the  easiest  way  to  measure  the  size  of  a  bubble  with  such  an  instrument  is  to 
make  use  of  near  critical-angle  scattering.  This  is  because  the  scattering  pattern  of  light 
from  a  bubble  in  water  contains  broadly  spaced  structures  which  are  indicative  of  the 
bubble  radius  a.  These  structures  are  present  as  the  scattering  angle  0  (relative  to  the 
direction  of  the  incident  beam)  decreases  below  the  critical  value  0C  =  82.8°.  The  original 

discovery  by  Marston  of  these  structures  is  described  in  Ref.  9  while  a  physical-optics 
approximation  (POA)  which  facilitates  a  simple  understanding  is  given  in  Ref.  1-5, 9. 

The  emphasis  of  the  research  supported  by  the  present  contract  (N00014-86-K- 
0242)  was  as  follows: 

1 .  Computations  that  showed  how  to  minimize  the  effects  that  natural  surface 
films  (on  bubbles  in  the  ocean)  will  have  on  scattering  data- the  procedure 
recommended  is  to  use  near-critical-angle  scattering  data  with  polarized 
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incident  light  with  the  detectors  in  the  plane  of  polarization.10’13  These 
calculations  indicate  that  there  may  be  an  optical  method  for  measuring  the  film 
thickness. 

2.  A  method  for  inversion  of  bubble  scattering  data14*15  on  an  event-bv-event 
basis— the  method  is  based  on  placing  detectors  at  the  critical  scattering  angle 
(0C  as  82.8  deg)  and  at  several  nearby  angles  0  <  0C.  The  bubbles  are 

assumed  to  be  illuminated  by  a  Gaussian  laser  beam. 

3.  Advice  given  on  technical  aspects  of  the  design  and  and  testing  of  an  optical 
bubble  spectrometer.15*16 

4.  Preliminary  consideration14  of  the  statistical  properties  of  the  scattering.  This 
may  yield  a  method  for  inversion  of  critical-angle  scattering  data  to  obtain  the 
size  spectrum  of  a  sample  of  events  (in  contrast  to  an  event-by-event 
inversion). 

In  the  present  Final  Report  some  aspects  of  items  1, 2,  and  4  will  be  noted  with  an 
emphasis  on  those  results  not  included  in  the  Technical  Report  (Ref.  10)  previously  issued 
and  available  from  D.T.I.C. 

In  the  Technical  Report10  different  values  of  the  refractive  index  for  ocean  water 
had  been  calculated  to  have  only  a  weak  effect  on  the  relevant  scattering  properties  of  the 
bubbles.  Consequently  examples  displayed  in  the  present  report  are  for  water  with  a 
refractive  index  of  4/3  »  1.3333. 

III.  RESEARCH  RESULTS  AND  ACCOMPLISHMENTS 

A.  Computational  Analysis  of  the  Effects  of  Surface  Films  on  the 
Optical  Scattering  Properties  of  Bubbles  in  Water 

Many  of  the  unresolved  issues  concerning  microbubbles  (of  natural  origin)  in  the 


ocean  concern  the  effects  of  adsorbed  films.  Such  films  are  thought  to  coat  most  (if  not  all) 
bubbles  in  sea  water  (see,  e.g..  Ref.  17).  The  films  should  affect  the  evolution  of  bubbles 


in  the  sea.  Recently  published  Soviet  data,18  for  example,  show  that  the  rate  at  which  a 
bubble  shrinks  in  sea  water  (when  subjected  to  pressure)  is  significantly  smaller  than  for 
bubbles  in  distilled  or  tap  water.  The  difference  is  attributed  to  the  supposed  presence  of 
surface  films  of  organic  substances.  The  film  "coats"  the  surface  of  the  gas  bubble. 
Glazman  suggests  (Sec.  I  of  Ref.  19)  that  the  film  thickness  may  range  from  10  nm  to  1 
pm. 

The  emphasis  of  the  present  research  was  to  examine  the  influence  such  films 
would  have  on  the  near  critical-angle  light  scattering  pattern  for  a  bubble  in  water.  The 
issue  is  important  since,  if  these  films  were  predicted  to  mask  the  critical-angle  pattern  for 
uncoated  bubbles,1'4*9  it  would  be  necessary  to  consider  major  design  changes  for  the 
proposed  bubble-size  spectrometer.  Fortunately,  it  was  found  that  the  films  should  have 
only  a  minor  influence  on  the  scattering  pattern  provided  certain  precautions  are  taken.10 
Of  these,  the  most  important  are:  (1)  to  make  use  of  the  coarse  structure  in  the  scattering 
pattern  but  to  average  over  the  fine  structure;  and  (2)  to  use  monochromatic  incident  light 
which  is  polarized  with  the  electric  field  parallel  to  the  scattering  plane,  that  is,  the  plane 
which  contains  the  detectors. 

The  method  used  in  the  computations  was  to  evaluate  the  exact  partial-wave  series 
solution  for  the  scattering  of  light  from  concentric  spheres  with  parameters  appropriate  for  a 
coated  bubble  in  water.  This  exact  solution,  given  by  Aden  and  Kerker  for  use  in  other 
scattering  problems20  will  be  referred  to  as  the  Aden-Kerker  Series  (AKS).  One  of  the 
important  accomplishments  of  the  present  research  was  to  develop  a  reliable  computer 
program  for  the  scattered  irradiance  as  a  function  of  scattering  angle  9  for  the  case  in 
which  the  sphere's  inner  radius  a  greatly  exceeded  the  wavelength  of  light  while  the 
coating  thickness  h  is  allowed  to  be  much  smaller  than  wavelength.  (It  appears  that  such  a 
reliable  program  was  not  previously  available.)  The  design,  testing,  and  a  listing,  of  the 
computer  program  are  summarized  in  Ref.  10. 


The  physical  situation  considered  is  shown  in  Fig.  1.  Light  is  incident  on  a  gas 
filled  bubble  coated  by  a  film  of  thickness  h  =  b  -  a  and  refractive  index  nc.  The 
refractive  index  of  the  surrounding  liquid  (water)  is  nw  =  4/3  while  the  refractive  index  of 
the  gas  within  the  bubble  is  n^  =  1.  Figure  1  also  shows  various  rays  reflected  from, 
and/or  refracted  by,  the  bubble.  Of  course  the  exact  AKS  is  based  directly  on  the  formal 
solution  of  Maxwell's  equations;  it  does  not  rely  on  geometrical  optics.  Nevertheless,  as 
will  be  discussed  below,  the  ray  diagram  is  a  first  step  toward  understanding  the  scattering 
pattern.  For  comparison,  the  dashed  rays  shown  are  drawn  so  as  to  illustrate  the  case  of  no 
coating  (which  is  equivalent  to  =  nw)  while  the  refraction  of  the  solid  rays  (at  the  water- 

coating  interface)  has  been  slightly  exaggerated  from  that  of  the  representative  case 

n<.  =  1.5.  This  was  done  so  as  to  show  the  direction  of  the  deflection  of  the  rays  introduced 

by  the  typical  coating. 

For  the  purposes  of  evaluating  the  AKS  partial-wave  series,  the  coating  index  nc 

was  usually  taken  to  be  1.5  since  that  is  close  to  the  refractive  index  of  many  oily  organic 
substances.  Calculations10  for  the  case  =  1.45  and  1.55  were  qualitatively  similar  to 
those  for  n^.  =  1.5  and  the  quantitative  differences  may  be  explained  in  terms  of  the  simple 
approximation  given  in  the  next  section.  The  scattered  irradiances  were  normalized  to  the 
irradiance  reflected  from  a  perfectly  reflecting  sphere  having  the  same  radius  a  as  that  of 
the  gas  pocket  within  the  bubble.10 

Figure  2  shows  the  normalized  irradiance  for  the  case  E-field  parallel  to  the 
scattering  plane.  The  bubble  radius  is  held  constant  at  a  =  7.5  pm  which  corresponds  to  a 
size  parameter  ka  =  100  (see  Ref.  10).  The  parameter  AOB  mentioned  in  the  caption  is  the 
ratio  a/b.  Figure  2(a)  shows  the  result  for  an  uncoated  bubble;  the  coarse  structure  (which 
gives  a  minimum  in  this  case  at  60  deg.)  is  clearly  seen  for  0  <  0C  =  82.8  deg.  Figures 

2(b)  -  (e)  show  the  increase  in  the  fine  structure  amplitude  which  occurs  as  the  coating 
thickness  h  =  b  -  a  increases.  Even  for  the  relatively  thick  coatings,  Fig.  2(f)  -  (i),  the 
coarse  structure  is  clearly  evident  and  is  similar  in  appearance  to  that  in  Fig.  2(a). 
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Figure  1 .  Coated  sphere  geometric  model.  Rays  are  specified  by  (pj)  where  p  ■  the 
number  of  chords  within  the  bubble,  and  /  ■  the  number  of  crossings  of  the  optic  axis. 
Dashed  lines  indicate  the  uncoated  sphere  ray  paths  while  solid  lines  represent  coated 
sphere  ray  paths.  Note  that  bending  of  the  rays  is  exaggerated  for  effect  The  subscript  b 
indicates  rays  which  pertain  to  the  outer  sphere  only. 
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Figure  2.  Evolution  of  the  fine  structure  as  a  function  of  coating  thickness  for  parallel 
polarized  imdiance  forks  *  100.  Coating  thickness  values  for  *  632.8  nm  are 

(a)  0  pin,  (b)  0.08  pm,  (c)  0.15  pm,  (d)  0.19  pm,  (e)  0.23  pm,  (f)  0.48  pm, 

(g)  0.66  pm,  (h)  0.93  pm,  and  (i)  1.13  pm.  The  corresponding  values  of  AOB  specified 
in  the  computation  were:  (a)  1.0,  (b)  0.990,  (c)  0.980,  (d)  0.975,  (e)  0.970,  (f)  0.940, 
(g)  0.920,  <h)  0.890,  (i)  0.870.  The  scattering  angle  runs  from  30  to  90  deg. 
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Figure  3  is  like  the  case  shown  in  Figure  2  except  that  the  full  angle  range  of  0  to 
180  deg  is  shown.  The  coarse  structure  is  clearly  broader  in  width  (and  easier  to  detect) 
than  the  relatively  narrow  forward  diffraction  lobes.  Indeed,  if  instead  of  a  =  7.5  jim  we 
take  a  =  75.5  pm,  the  angular  quasi-period  of  the  forward  diffraction  lobes  decreases  by  a 
factor  of  10  while  the  coarse  quasi-period  decreases  by  a  factor  =  3;  consequently,  the 
difference  in  quasiperiods  is  even  more  clearly  seen  for  larger  bubbles.  Inspection  of  Fig. 
3(b)  shows  that  for  this  case  where  h/a  has  a  relatively  large  value  (1. 1/7.5  =  0.15)  the 
coarse  structure  manifests  a  noticeable  shift  towards  larger  scattering  angles. 

Figure  4  shows  the  scattered  irradiance  for  polarized  light  also  for  a  case  in  which 
the  E  field  is  parallel  to  (i.e.,  lies  within)  the  plane  of  scattering.  In  all  of  the  curves  a  = 
37.8  pm.  The  solid  curve  is  for  a  coated  bubble  with  h  =  1  pm  and  was  given  by  the 
exact  AKS  partial-wave  series.  The  curve  with  short  dashes  is  for  an  uncoated  bubble  and 
was  given  by  the  Mie  series.  These  numerical  values  for  h  and  a  are  for  the  case  of  red 
light  having  a  wavelength  in  air  of  633  nm.  The  slowly  varying  curve  with  intermediate 
length  dashes  is  the  result  of  the  previous  physical-optics  approximation  (POA)  for 
scattering  from  a  homogeneous  bubble  previously  formulated  by  Marston  and 
Kingsbury.1*2  That  curve  is  seen  to  approximate  the  coarse  structure  in  the  Mie  curve  (as 
well  as  the  data  for  scattering  from  uncoated  bubbles,  see  Ref.  5).  The  right-most  smooth 
curve  (which  alternates  a  long  dash  with  a  short  dash)  is  a  shifted  plot  of  the  POA  curve. 
That  shift  was  calculated  geometrically  by  Marston  so  as  to  obtain  an  approximate 
description  of  the  scattering  from  a  coated  bubble  near  the  critical  scattering  angle. 

Figure  4  shows  that  the  principal  effect  of  the  coating  on  the  coarse  structure  is  to 
shift  slightly  its  angular  location  in  the  direction  of  increased  scattering  angle.  Notice, 
however,  the  coarse  structure  should  not  be  masked  by  the  presence  of  the  coating.  Hence 
size  spectrometers  which  make  use  of  the  coarse  structure  should  still  function  if  the  bubble 
is  coated.  The  coating  may,  however,  introduce  a  small  error  (or  shift)  in  the  bubble  size 
estimated  from  scattering  data.  It  is  noteworthy  that  the  effects  of  the  coating  on  the 
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Figure  3.  Calculated  normalized  scattering  irradiances  I2  for  ka  =  100,  nw  =  4/3,  = 

1.5, 9  =  0°  to  180°  for  the  log  of  the  parallel  polarization.  The  solid  curve  is  from  Mie 

theory.  The  thin-dashed  curve  is  the  calculated  coated  sphere  result  for  (a)  h  =  0.26  fim, 
(b)  h  =  1.14 11m.  All  values  of  h  are  for  HeNe  laser  light  in  water  with  =  632.8  nm. 

The  inner  radius  a  =  7.5 11m. 
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orthogonal  polarization  (E-field  perpendicular  to  the  scattering  lane)  were  much  more 
significant;  that  choice  of  polarization  may  be  unsuitable  for  use  in  a  size  spectrometer. 
These  polarization  effects  may  be  explained  in  terms  of  physical  arguments.10 

Figure  5  is  for  the  same  radius  and  coating  thickness  as  Fig.  4  but  it  shows  the 
normalized  irradiance  for  the  case  in  which  the  E-field  is  perpendicular  to  the  plane  of 
scattering.  The  dashed  curve  is  the  case  of  an  uncoated  bubble  and  is  given  by  Mie  theory 
while  the  solid  curve  is  the  coated  bubble  and  is  given  by  the  AKS.  The  effect  of  the 
coating  on  the  coarse  structure  is  more  than  that  of  the  shift  evident  in  Fig.  4. 

The  Technical  Report  (Ref.  10)  presents  computed  irradiance  plots  for  bubbles 
having  radii  as  large  as  1 89  (J.m.  The  effects  of  optical  absorption  by  the  coating  material 
were  also  explored.  It  turns  out  that  unless  the  coating  contains  a  dye,  absorption  should 
not  have  a  major  effect  This  is  because  of  the  small  thickness  (h  £  1  pm)  of  the  anticipated 
films.  As  noted  in  the  supplement  at  the  end  of  that  report,  some  of  POA  (physical  optics 
approximation)  plots  presented  were  subject  to  small  computational  errors.  These  errors 
were  subsequently  corrected  and  are  not  present  in  the  POA  plot  shown  here  in  Fig.  4. 

B.  Physical  Optics  Approximation  for  the  Scattering  from  a  Coated 
Bubble-Explanation  of  the  Shift  in  the  Coarse  Structure  and 
Polarization  Effects 

Figures  3  and  4  show  that  the  coating  shifts  the  coarse  structure  to  larger  scattering 
angles.  To  understand  this  shift,  it  is  appropriate  to  recall  the  physical  origin  of  the  coarse 
structure  in  the  physical-optics  approximation  for  the  case  of  an  uncoated  bubble.1’2’9  The 
structure  arises  because  of:  (a)  the  effects  of  diffraction  on  the  wavefront  reflected  in  the 
vicinity  of  the  critical  angle  of  incidence,  and  (b)  the  interference  pattern  of  the  wavefront 
associated  with  transmitted  rays  (predominantly  the  one-chord  ray)  with  the  reflected 
wavefront.  Diffraction  is  particularly  significant  when  calculating  the  scattering  amplitude 
associated  with  the  reflected  wavefront  (a)  because  of  the  abrupt  drop-off  of  the  local 


scattering  angle  (deg) 


reflection  coefficient  as  the  angle  of  incidence  i  falls  below  the  critical  value  = 
sinfnw’1).  For  the  uncoated  bubble  this  reduction  occurs  when  the  scattering  angle  6  is 


18 


b: 


p 


increased  above  the  critical  value  0C 

0  =  7t-2i  =  Jt  -  sin'^nj )  s  82.8°  (1) 

c  c  w 

The  diffraction  and  interference  patterns  caused  by  effects  (a)  and  (b)  are  present  in  the 
region  0  <  0C.  These  patterns  have  about  the  same  quasi-period  and  they  manifest 

themselves  as  the  coarse  structure  in  Fig.  2(a)  and  in  the  Mie  curve  in  Fig.  4.  The 
interference  of  the  aforementioned  ray  amplitudes  with  that  of  the  far- side  ray,  labeled  (2,1) 
in  Fig.  1,  gives  rise  to  the  fine  structure.  We  will  limit  our  attention  to  the  coarse  structure 
since  it  is  to  be  used  to  determine  the  size  of  the  bubble. 

Now  what  happens  when  you  coat  the  bubble?  The  relevant  interface  at  which 
there  is  an  abrupt  transition  to  total  reflection  becomes  the  inner  surface  of  the  coating. 
Denote  the  local  angle  of  incidence  for  a  ray  at  this  surface  by  ia  where  the  subscript  a 

denotes  the  radius  of  the  relevant  interface.  There  will  be  an  abrupt  transition  to  total 
reflection  here  when  ia  £  i^  where  the  new  critical  angle  is 

‘ca  =  (2) 

Now  the  ray  which  reflects  from  the  inner  surface  will  be  refracted  by  the  coating-water 
interface  and  leaves  the  bubble  with  a  well  defined  scattering  angle  relative  to  the  direction 
of  the  incident  beam.  For  that  ray  having  ia  =  ica,  denote  the  resulting  scattering  angle  by 
0C'.  Evidently  0C'  is  the  critical  scattering  angle  of  the  coated  bubble.  The  presence  of  the 
coating  shifts  the  critical  scattering  angle  by  an  amount 

eJ-ec  =  A,  O) 

Now  A  is  a  function  of  the  coating-thickness  to  radius  ratio  h/a  as  well  as  nw  and  nc. 
Exact  calculation  of  A  requires  the  solution  of  a  set  of  transcendental  equations;  however, 


Marston  obtained  the  leading  term  in  a  power  series  expansion  of  A(h/a).  The  resulting 
approximation  is 


A  = 


-  (n|  -  l)1/23  +  0(h/a)2, 


(4) 


where  0(h/a)2  is  a  function  which  vanishes  as  rapidly  as  (h/a)2  as  h  — » O.  The 
derivation  of  Eq.  (4)  was  repeated  by  Cleon  Dean,  a  graduate  student  partially  supported  by 
this  contract.  Notice  that  A  -4  O  if  the  coating  index  ->  nw  or  if  the  thickness 


h  -4  0. 

To  approximate  the  effects  of  the  coating  on  the  coarse  structure,  the  POA  for  an 
uncoated  bubble  may  be  evaluated  with  the  angular  argument  shifted  by  A.  Let 
IJPOA(0,ka)  be  the  normalized  irradiance  at  scattering  angle  9  as  calculated  (for 

polarization  j)  by  the  original  POA  for  an  uncoated  bubble.1’2’5  The  resulting 
approximation  for  the  scattering  from  a  coated  bubble  is 


tcpoa.qi  .  .POA.f.  .  .  . 
L  (0,ka)  »  Ij  (0  -  A,  ka). 


(5) 


where  the  coating  parameters  n<.  and  h/a  enter  into  this  approximation  only  through  the 
approximation  of  A  by  the  first  term  in  Eq.  (4).  The  main  idea  behind  Eq.  (5)  is  that  the 
principal  effect  of  a  thin  coating  (at  least  for  polarization  j  =  2  as  noted  below)  is  to  shift 
the  scattering  angles  of  the  critical  ray  and  the  single-chord  ray  in  Fig.  1  by  similar 
amounts.  This  shift  results  in  a  corresponding  shift  of  the  diffraction  and  interference 
pattern. 

How  well  does  Eq.  (5)  do  in  approximating  the  exact  AKS  results?  Consider  first 
polarization  j  =  2  which  is  the  case  of  E  field  parallel  to  the  scattering  plane.  The  curve 
in  Fig.  4  with  alternate  long  and  short  dashes  gives  l2CPOA(6,  ka  =  500).  It  is  evident  that 
this  CPOA  describes  the  general  rise  in  irradiance  of  the  first  coarse  peak  as  9  decreases 
from  85  to  75  deg.  The  shift  in  the  location  of  the  peak  near  65  deg  is  also  well 
approximated  by  the  model.  This,  and  several  other  comparisons  done  for  j  =  2  with  ka 
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from  100  to  1000  show  that  the  CPOA  does  about  as  well  in  approximating  the  coarse 
structure  for  a  coated  bubble  as  the  POA  does  for  an  uncoated  bubble.  See  Fig.  6. 

Now  what  about  polarization  j  =  1  which  is  the  case  of  E-field  perpendicular  to  the 
scattering  plane?  It  was  evident  from  Fig.  5  (and  other  plots1®)  that  the  effect  of  the  coating 
on  the  coarse  structure  is  more  than  just  an  angular  shift.  The  modulations  of  the  coarse 
structure  introduced  by  the  coating  may  be  attributed  to  the  reflection  of  light  at  the  coating 
water  interface)  as  illustrated;  e.g.,  by  ray  (0,0)^  in  Fig.  1.  The  scattering  amplitude 

associated  with  this  ray  is  larger  for  j  =  1  than  for  j  =  2  since  for  polarization  j  =  2  the 
local  angle  of  incidence  (ib  in  Fig.  1)  lies  close  to  the  Brewster  angle  of  incidence  iB 
=  tan'1(nc/nw)  for  this  interface.  For  a  plane  surface,  the  Fresnel  reflection  coefficient 
vanishes  when  ib  -  iB  for  polarization  j  =  2  but  for  j  =  1  there  is  no  minimum  in  the 
reflection  coefficient.  Hence  we  understand  why  for  j  =  1  the  coarse  structure  is  disrupted 
by  the  coating  so  that  scattering  with  j  =  I  may  be  difficult  to  determine  the  bubble’s 
radius. 

Thus  far  we  have  discussed  only  the  usual  case  of  nc  >  nw  =  1.333  where  the 
specific  examples  in  Fig.  4  and  5  had  nc  =  1 .5.  Inspection  of  Eq.  (4)  leads  to  the 
prediction  that  A  <  0  in  the  unusual  case  of  ^  <  nw.  To  confirm  our  general 

understanding  of  the  effects  of  coatings,  we  computed  the  exact  AKS  for  the  case  j  =  2,  a 
=  37.5  pm,  h  =  3.3  pm,  and  nc  =  1.25.  Here  the  coarse  peak  was  shifted  to  smaller  0 

(i.e.,  A  <  0  as  predicted)  and  Eq.  (5)  well  approximated  the  general  rise  in  the  coarse  peak 
as  0  decreased  below  85  deg.  It  should  be  noted,  however,  that  bubbles  in  the  ocean 
would  not  have  coating  for  which  nc  <  nw. 

C.  Brewster- Angle  Scattering  from  Coated  and  Uncoated  Bubbles: 
Applications  to  Measurement  of  Film  Thickness  and  to  the 
Discrimination  of  Bubbles  from  Particles 
The  present  section  is  concerned  with  properties  of  scattering  of  polarized  light 
having  the  E-field  parallel  to  the  scattering  plane  (j  =  2).  Inspection  of  the  Mie  theory 
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PLOTS  AS  IN  FIG.  4  EXCEPT  THAT  THE  COATED  BUBBLE  MODELS  HAVE 


FIGURE 


results  for  the  uncoated  bubbles  displayed  as  the  solid  curves  in  Fig.  3(a,b)  reveals  the 
following  major  features  when  the  scattering  angle  0  =  106  deg: 

(a)  Let  <I2>  denote  a  local  angular  average  of  the  irradiance  which  masks  the  fine 
structure.  This  may  be  called  the  baseline  irradiance.  It  is  evident  the  <l2>  is 
minimized  for  0  close  to  106  deg. 

(b)  The  visibility  of  the  fine  structure  fringes  in  the  vicinity  of  a  given  scattering 
angle  0  may  be  denoted  by 


Vis(0)  = 


^max  ’  *min 

I  +T“  ’ 

max  min 


(6) 


where  Imax  and  1,^  are  local  maximum  and  minimum  values  of  I2  closest  to  the 
specified  0.  Inspection  of  the  Mie  results  (solid  curves  in  Fig.  3(a,b))  shows  that  Vis(0) 
has  a  local  minimum  for  0  close  to  106  deg. 

Consider  the  local  angle  of  incidence  i  for  the  ray  reflected  from  the  water-air 
interface  of  the  uncoated  bubble.  The  aforementioned  features  are  a  consequence  of  i 
being  close  to  the  relevant  Brewster  angle  iB  =  tan^fiv1)  =  tan-^OJS)  =  36.9  deg.  The 
scattering  angle  corresponding  to  i  =  iB  is 


0R  =  180  deg  -  2ig  =  106.2  deg,  (7) 

which  may  be  called  the  Brewster  scattering  angle.  Since  the  plane  surface  reflection 
coefficient  has  a  null  when  i  =  iB,  it  is  to  expected  that  <l2>  is  minimized  for  0  =  0B  as 

was  previously  noted  in  Ref.  3  and  21. 

Now  the  fine  structure  fringes  in  the  angular  region  between  0  of  82.8  deg  and 
120  deg  are  primarily  a  consequence  of  the  interference  between  the  far-side  ray  (the 
dashed  (2,1)  ray  in  Fig.  1)  with  the  reflected  ray  (the  dashed  (0,0)  ray).  Since  the 
amplitude  of  the  reflected  ray  is  minimized  near  0  =  0B,  so  is  the  visibility,  Eq.  (6),  of  the 

interference  pattern. 


Figure  3  shows  changes  in  the  scattering  near  106  deg  which  result  from  coating 
the  bubble  with  a  film  having  an  index  ^  =  1.5.  These  results  are  representative  of  other 

examples  plotted  in  Ref.  10.  The  principal  change  is  an  increase  of  the  baseline  irradiance 
<l2>  near  lOo  deg.  This  change  is  to  be  expected  since  the  coating  will  mask  the  local  null 
in  the  Fresnel  reflection  coefficient.  Notice  that  <12(106  deg)>  is  greater  for  the  thicker  of 

the  two  coatings  (Fig.  3(b)  as  opposed  to  Fig.  3(a)).  Results  displayed  in  Ref.  10  suggest 
that  <l2(106deg)>  is  a  monotone  increasing  function  of  the  film  thickness  h  for  a  range 

of  h;  however,  due  to  limited  resources,  this  problem  has  not  been  examined  in  detail.  It  is 
plausible  that  the  film  thickness  h  could  be  determined  from  measurements  of 
<l2(106deg)>  or  from  the  ratio  p  discussed  below. 

Another  change  evident  in  Fig.  3  caused  by  the  coating  is  in  the  amplitude  of  the 
fine  structure  near  106  deg.  These  changes  would  certainly  affect  the  visibility  function, 

Eq.  (6)  near  106  deg  though  they  might  not  be  as  useful  for  determining  the  film  thickness. 

There  is  an  additional  application  of  this  Brewster-angle  scattering  which  should  be 
noted.  Suppose  bubbles  in  the  ocean  were  not  coated  or  that  as  in  Fig.  3(a),  the  coating 
was  so  thin  as  to  have  only  a  minor  effect  on  <12  (106  deg)>.  Then  it  would  be  an  easy 

matter  to  verify  that  the  scatterer  was  in  fact  a  bubble  by  measuring  the  ratio 

p  =  <12(82.8  deg)>/<l2(  106  deg)>.  (8) 

For  such  a  clean  bubble,  p  will  be  nearly  independent  of  the  radius  a  of  the  gas  pocket  if 
a  >  7  pm.  Computational  results  suggest  that  p  will  lie  between  25  and  100.  If  the 
scatterer  was  a  small  spherical  particle  or  biological  cell  having  a  refractive  index  >  nw,  Mie 
theory  typically  gives  p  <  1.  Hence  such  non-bubble  scattering  events  could  clearly  be 
identified.  Now  for  coated  bubbles,  it  appears  that  p  decreases  and  can  have  values  as 
small  as  p  =  5.  Hence  the  identification  of  non-bubble  events  from  the  measured  value  of 


D.  Detector  Placement  and  Inversion  of  Scattering  Data  to  Obtain  the 


Bubble  Size  on  an  Event-by-Event  Basis 
In  addition  to  the  task  of  computing  and  understanding  the  optical  effects  of  surface 
films,  the  other  major  task  of  the  contract  was  to  consider  how  the  light  detectors  should  be 
placed  so  as  to  facilitate  a  measurement  of  a  bubble's  radius  a.  The  geometry  considered 
was  intended  to  be  realizable  on  instruments  of  the  Wyatt  Technology22  design  originally 
intended  for  studying  the  scattering  from  single  marine  microparticles  in  situ.23  A  region 
of  water  is  considered  to  be  illuminated  by  a  Gaussian  beam  of  polarized  monochromatic 
light.  For  such  a  beam,  the  incident  irradiance  has  the  spatial  dependence 


‘inc  =  *o  exP<-rW). 


(9) 


where  b  is  a  parameter  which  specifies  the  effective  beam  radius  and  r  is  the  distance 
from  the  axis  of  the  beam.  Bubbles  drift  through  the  scattering  region  and  the  light 
scattered  by  the  bubbles  is  collected  by  a  set  of  detectors;  each  having  a  well  defined  central 
scattering  angle  9  relative  to  the  direction  of  the  incident  beam.  The  basic  problem  is  to 
consider  how  to  select  the  9  of  a  set  of  detectors. 

As  described  in  Sec.  ELA-C  of  this  report,  the  effects  of  a  thin  coating  on  the 
scattering  pattern  of  a  bubble  are  expected  to  be  small  if  the  incident  light  is  polarized  with 
the  E-field  parallel  to  the  plane  of  the  detectors.  Consequently,  we  will  limit  our  attention 
to  such  a  configuration  and  model  the  scattering  as  if  the  bubble  were  not  coated.  In  an 
experiment,  the  error  introduced  by  a  coating  in  the  data  inversion  process  could  be 
estimated;  however,  it  is  presently  anticipated  that  coating  would  not  be  the  principal  source 
of  error  when  determining  either  the  bubble  radius  or  the  absolute  number  density  of 
bubbles. 

Each  detector  produces  a  voltage  proportional  to  the  scattered  irradiance  averaged 
over  the  range  of  angles  viewed.  Denote  this  voltage  by  vj(9,t)  where  9  is  the  central 

scattering  angle  of  that  detector  and  the  polarization  index  j  has  the  value  2  for  the 


configuration  considered.  The  vj  depend  on  time  t  because  the  bubble  moves  through 
the  beam.  The  Vj(0,t)  are  assumed  to  be  recorded  as  a  function  of  time.  The  principal 
problem  with  the  use  of  a  Gaussian  beam,  Eq.  (9),  is  that  the  experimenter  does  not  know 
how  close  the  bubble  comes  to  the  axis  of  the  beam.  Let  the  minimum  distance  of  the 
center  of  a  bubble  from  this  axis  be  denoted  by  y.  Define  an  "event"  as  the  passage  of  a 
single  bubble  through  the  beam;  y  may  be  thought  of  as  an  impact  parameter  relative  to  the 
beam's  axis  for  a  given  event.  If  the  radius  a  of  the  bubble  is  «  b,  the  detected  voltage 
will  be  maximized  if  y  =  0  (such  that  the  bubble  is  centered  on  the  beam).  Now  for  such 
an  event,  the  detector  at  the  critical  angle  0  =  0C  =  82.8  deg  will  produce  a  voltage  V2 

which  increases  monotonicallv  with  the  bubble  radius  a.  It  would  then  be  a  simple  matter 
to  invert  the  data  on  V2  to  obtain  the  radius  for  the  bubble  in  each  event  Unfortunately, 
when  y  is  not  known,  inversion  of  the  V2  data  is  complicated. 

Nevertheless,  a  procedure  was  formulated  whereby  the  bubble  radius  for  each  event 
could  be  determined  within  a  narrow  range.  It  is  assumed  that  a«b  and  that  the 
instantaneous  scattering  pattern  could  be  well  approximated  by  that  for  a  plane  wave  with 
an  incident  irradiance  given  by  Eq.  (9)  with  r  taken  to  be  the  instantaneous  distance  of  the 
bubble  from  the  beam  axis.  Some  noteworthy  features  of  this  procedure  are: 

(a)  Detectors  are  to  be  placed  at  the  critical  scattering  angle,  0C  =  82.8  deg  and  at 
several  nearby  angles  having  0  <  0C. 

(b)  The  voltages  wj  fr°m  these  detectors  will  occur  in  ratios  that  are  somewhat 

slowly  varying  functions  of  bubble  radius  a  because  the  coarse  structure 
(which  occurs  in  the  region  0  <  0C)  varies  slowly  with  a.  Specifically,  the 
ratios  V2(0,t)/v2(0  =  0c,t)  should  be  considered  for  that  time  t  which 
maximizes  V2(0c,t)  for  a  given  event. 


(c)  It  is  anticipated  that  only  a  narrow  range  of  radius  values  will  give  theoretical 
values  for  the  aforementioned  ratios  consistent  with  data  for  scattering  from 
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single  bubbles.  The  mean  value  in  this  range,  call  it  a^p  may  be  taken  as  the 


measured  value  for  that  event 

(d)  The  inversion  procedure  may  be  subjected  to  various  statistical  and  internal 
consistency  checks  to  test  if  it  is  operating  properly.  For  example,  from  a^p 
and  the  V2(0c,t)  data  one  could  work  backwards  to  determine  the  impact 
parameter  y  of  each  event.  The  resulting  y  should  be  randomly  distributed 
in  the  viewed  region.  (Note,  the  viewed  region  of  water  may  be  limited  by 
cutting  off  the  irradiance  so  that  iinc  drops  to  zero  abruptly  for  large  values  of 


The  details  of  this  procedure  were  given.14  Rather  than  repeat  the  description,  the 
relevant  portion  of  that  report  is  reproduced  as  Appendix  A  to  the  present  Final  Report 
Before  proceeding,  the  reader  is  asked  to  examine  Appendix  A.  Table  II  of  the  appendix 
lists  recommended  detector  locations  and  angular  widths. 

Some  comments  on  the  method  of  event-by-event  data  inversion  presented  in 
Section  A  of  Appendix  A  are  appropriate. 

(a)  Figures  A7-A10  of  A,ppendix  A  were  claimed  to  be  plots  of  the  normalized 
irradiances  as  a  function  of  the  bubble  size  parameter  ka  (see  Table  I  for 
conversion  to  radius  values)  as  given  by  the  physical-optics  approximation. 
Since  those  plots  were  prepared,  a  minor  error  in  the  POA  computer  algorithm 
used  was  corrected  (see  the  Supplement  Section  of  Ref.  10).  The  corrections 
are  minor  as  may  be  seen  by  comparing  the  corrected  POA  curve,  Fig.  7,  with 
Fig.  A8.  These  curves  given  the  normalized  scattering  irradiance  I2  for  0  = 

75  deg.  The  corrections  do  not  affect  the  applicability  of  the  procedure.  (Note: 
the  original  presentation  and  plots  of  the  POA,  Ref.  1-3  were  free  from  the 
aforementioned  computational  error  as  are  Fig.  4, 6,  A2  of  the  present  report.) 

(b)  Figures  A7-A10  were  intended  to  illustrate  the  method  of  data  inversion.  If 
this  method  is  actually  used  for  a  real  instrument,  it  would  probably  be 


1 


IS  PROPORTIONAL  TO  I0(ka) 
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necessary  to  prepare  plots  (or  tables)  of  the  voltage  ratios  V2(0)/v2(8  =  0C)  for 
the  actual  detector  widths  used  by  averaging  the  Mie  theory  results  over  the 
appropriate  angular  region  for  the  range  of  bubble  radii  encountered  in  the 
experiments. 

(c)  To  determine  the  absolute  bubble  concentration  it  is  necessary  to  know  the 
cross  sectional  area  of  the  scattering  region  through  which  the  water  flows. 

The  significance  and  control  of  this  area  is  discussed  in  Sec.  B.3  of  Appendix 
A. 

All  of  the  aforementioned  analyses  assumes  that  there  is  only  a  single  bubble  in  the 
scattering  volume  during  the  time  t  when  Vj(0c,t)  has  its  maximum  value.  It  is  also 

assumed  there  is  no  significant  scattering  of  light  bubbles  in  the  volume  of  water  which  the 
beam  passes  through  before  it  reaches  the  scattering  volume.  The  problem  of  how  to  avoid 
or  minimize  the  effects  of  multiple  scattering  was  addressed^  and  the  suggestions  will  not 
be  repeated  here.  It  is  noteworthy  that  data  which  does  not  fit  the  predicted  pattern  for 
single  bubbles,  presented  here  in  Appendix  A  and  in  Sec.  C  (for  the  Brewster  scattering 
angle  data),  could  be  attributed  either  to  scattering  from  more  than  one  bubble  or  from 
particles  and/or  biological  cells. 

E.  Statistical  Properties  of  the  Photodetector  Voltages  and  a 
Possible  Alternative  to  Event-by-Event  Data  Inversion 

The  maximum  in  the  detector  voltage  output  for  a  given  bubble  event  will  depend 
on  how  close  the  bubble  came  to  the  axis  of  the  beam;  this  minimum  distance  is  the  impact 
parameter  y  discussed  in  the  previous  section.  For  a  specific  photodetector,  the  voltage 
output  will  pass  through  a  maximum  value  as  the  bubble  passes  through  the  beam,  see  Fig. 
A1  of  the  Appendix.  Let  v  denote  this  maximum  value  for  a  given  event;  v  depends  on 
the  impacted  parameter  y  of  the  event  and  for  some  value  of  y,  v  takes  on  a  maximum 
value  v0  which  will  depend  on  the  radius  a  of  the  bubble.  Assumptions  will  be  noted 
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below  concerning  the  statistics  of  y  and  the  scattering  of  light  from  small  bubbles  in  a 
Gaussian  beam.  These  assumptions  facilitate  the  prediction,  below,  of  the  statistical 
properties  of  the  random  variable  z  =  v/v0  for  a  sample  of  events  in  which  all  bubbles 

have  the  same  radius.  For  the  usual  case  in  which  bubbles  in  a  sample  of  events  are  not  all 
of  the  same  size,  that  result  could  be  used  to  predict  the  probability  density  function  F(v) 
of  the  absolute  voltages;  it  is  also  necessary,  however,  to  know  the  size  distribution  (per 
radial  interval)  n(a)  of  the  polydisperse  population  of  bubbles  in  the  sample.  It  appears 
that  in  principle  the  problem  may  be  turned  around  so  that  an  unknown  n(a)  might  be 
found  from  the  observed  F(v);  however,  the  details  and  stability  of  such  an  inversion 
process  are  yet  to  be  examined. 

The  relevant  statistical  properties  of  z  =  v/v0  for  a  monodisperse  size  distribution 
of  bubbles  will  now  be  derived.  It  is  assumed  that  the  irradiance  profile  is  described  by 
Eq.  (9)  and  that  a  «  b  so  that  the  scattered  irradiance  for  a  plane  wave  of  irradiance  given 
by  the  local  value  at  the  bubble's  center, ^  these  assumptions  given 

v(y)  =  vQ  exp(-y2/b2),  (10) 

where  v0  =  v(y  =  0).  Now  the  impact  parameter  y  of  each  event  is  a  random  variable.  It 
will  be  assumed  that  y  is  uniformly  distributed  between  y  =  0  and  y  =  ymax  where  for 
now  we  take  ymax  — >  «*>.  Let  dN  denote  the  number  of  events  in  interval  between  y  and 
y  +  dy.  The  uniform  distribution  of  impact  parameters  gives 


(11) 


where  D  is  a  constant  (i.e.,  D  does  not  depend  on  y). 

The  random  variable  z  =  v/v0  is  related  to  y  byEq.  (10).  The  desired  quantity  is 

the  probability  density  function  F(z).  For  a  sample  of  events  having  a  large  number  of 
bubbles  of  the  same  size,  the  frequency  of  occurrence  of  event  voltages  between  v  and 
v  +  dv  will  be  proportional  to  F(v/v0)v0'Mv.  The  expression 
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F(z)  =  v0ldN/dvl  =  (dN/dy)ldy/dzl,  (12) 

may  be  shown  to  be  consistent  with  rules  governing  the  statistics  of  transformed  random 
variables.  Now  Eq.  (10)  gives 

dz/dy  =  -  2yb'2z  =  -2b  1  z[-ln(z)],  (13) 

where  0  £  z  <.  1.  Equations  (1 1)  -  (13)  yield 

F(z  ^  1)  =  A[z(-ln  z)1/2]  1  (14a) 

F(z  >  1)  =  0  (14b) 

where  the  constant  A  =  Db/2. 

The  density  function  F  in  Eq.  (14)  diverges  as  z  —»  0  because  the  Gaussian  beam 
was  assumed  to  be  unbounded  and  upper  bound  ymax  for  the  impact  parameters  was  also 

taken  to  be  infinite.  For  a  beam  in  which  the  irradiance  is  cut  off  at  some  finite  distance 
from  the  axis,  this  divergence  in  F(z)  is  removed.  As  z  — >  l,F(z)  diverges  because  of 
the  flatness  of  the  irradiance  profile  in  the  center  of  the  beam.  Except  for  these  divergences 
F  in  Eq.  (14)  is  smooth  and  F  takes  on  a  minimum  value  at  some  intermediate  z.  The 
divergence  of  F  as  z  — >  1  is  desirable  as  noted  below. 

Let  F(v)  denote  the  probability  density  function  of  the  absolute  voltage  v  fora 
sample  consisting  of  many  events  in  which  the  bubbles  were  not  all  of  the  same  size.  The 
frequency  of  occurrence  of  event  voltages  between  v  and  v  +  dv  will  be  proportional  to 
F(v)  dv.  Let  n(a)  denote  the  number  of  bubbles  having  radii  between  a  and  a  +  da;  n(a) 
is  the  size  distribution  for  the  sample.  The  statistical  properties  of  the  v  are  governed  by 

oo 

F(v)  =  ft  1  J  n(a)  F(v/vQ)da,  (15) 

o 


31 


where  N  is  the  total  number  of  bubbles  in  the  sample;  v0  depends  on  the  bubble  radius  a 
through  the  scattering  law  appropriate  for  the  detector  considered.  It  is  desirable  that  v0 
be  a  monotonic  increasing  function  of  a;  this  may  be  achieved  by  letting  the  detector  be 
centered  on  the  critical  scattering  angle  9  =  82.8  deg  and  by  polarizing  the  incident  beam 
so  that  the  E-field  lies  in  the  scattering  plane  (i.e.,  j  =  2).  For  example,  the  Physical  Optics 
Approximation1  •2*5  yields  the  following  proportionality 

v  =  Ga2  (16) 

o 

where  G  is  a  constant  which  depends  on  the  light  power  and  the  detector  sensitivity.25 
Notice  that  divergence  in  F(z)  as  z  — >  1  is  desirable  since  this  causes  F(v)  to  retain 
information  concerning  the  size  distribution  n(a).  Furthermore,  it  would  be  desirable  to 
remove  the  divergence  in  F(z)  as  z  — >  0  by  cutting  off  the  beam  at  some  finite  distance 
from  the  axis.  Of  course,  for  a  finite  beam  with  a  flat  irradiance  profile  F(z) «  5(z- 1) 
where  8  is  the  Dirac  8  function;  in  that  case  the  inversion  to  give  n(a)  from  data  for 
F(v),  by  way  of  Eqs.  (15)  and  (16),  would  be  easy. 
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VI.  APPENDIX  A:  DETAILS  OF  INVERSION  PROCEDURE  (REPRODUCED  FROM  REF.  14) 

A.  Detector  Placement  and  Proposed  Inversion  Procedure  Based  on  Relative  Scattered 
Irradiancet 


1.  Introduction— 

In  item  4(0  of  my  letter  of  April  25  I  indicated  that  other  detectors  would  be  needed  at 
angles  more  forward  than  6C  =  82.8  deg.  I  now  explain  the  positioning  and  use  of 
such  detectors.  Let  vj(0,t)  denote  the  voltage  output  of  a  given  detector  multiplied  by 
a  sensitivity  factor  determined  during  calibration:  j  =  2  for  a  detector  in  the  plane 
parallel  to  the  incident  wave's  E  polarization;  j  =  1  for  a  detector  in  the  plane 
perpendicular  to  the  incident  wave's  E  polarization;  8  denotes  the  central 
scattering  angle  viewed  bv  the  detector.  For  a  perfectly  isotropic  and  unpolarized 
point  source  located  in  the  center  of  the  scattering  volume,  the  vj  for  all  8  should 


have  the  same  value. 

Suppose  now  the  vj(8,t)  are  in  response  to  a  bubble  which  drifts  through  the 
scattering  volume.  As  a  function  of  time  t  the  signal  for  an  "event”  might  look 
something  like  that  shown  in  Fig.  1: 


Fig.  A1 


where  I  have  chosen  to  take  0  =  82.8°  and  j  =  2  since  that  is  a  particularly 
convenient  reference  for  reasons  noted  below.  The  Gaussian  shape  of  curve 
reflects  the  assumption  that  the  bubble  moves  at  a  fixed  velocity  through  a  Gaussian 
beam;  however,  neither  assumption  is  essential  to  the  analysis  which  follows.  Let 


VjCe.W  (82.8  °))  ■  Vj(9)  such  that,  for  now,  I  only  use  data  from  all  the  detectors  for 
the  time  interval  in  which  v2(82.8°,t)  i9  maximized,  that  is,  at  t  =  (82.8°). 

(Note  it  is  to  be  anticipated  that  for  spheres  the  other  vj(9,t)  would  be  maximized  at 

the  same  t  if  the  bubble  radius  a  « the  beam  width.)  Unfortunately,  the  values  of 
v2  (82.8°)  ■  v2  (82.8,  tnjg*  (82.8°))  do  not  by  themselves  give  the  bubble  radius  on  a 

event-by-event  basis  since  the  impact  parameter  y  of  a  bubble,  relative  to  the  center 
of  the  beam,  is  not  known.  Nevertheless.  I  show  in  Section  B  that  the  absolute 
magnitude  of  v2  (82.8°)  for  each  event  is  useful  both  for  augmenting  and 


checking  the  data  inversion  procedure  given  below. 

The  problem  now  reduces  to  finding  the  bubble  radius  a  for  a  given  event 
from  the  relative  outputs 

Wj(6)  -  Vj(9yv2(82.80)  (1) 

for  a  small  set  of  detectors  at  various  9.  The  procedure  listed  below  is  based  on 
ideas  which  I  presented  at  NORDA  in  January  1985. 

Renew  of  Scattering  Theory  fieaalia- 

Figures  2-5  show  the  computed  normalized  irradiances  I2  for  j  =  2  (E  field  parallel 
to  the  plane  of  scattering).  The  normalization  is  my  usual  one  in  which  the 
geometric  (isotropic)  reflection  from  a  perfectly  reflecting  sphere  of  the  same 
radius  gives  Ij  =  1.  Except  as  noted  below  the  solid  curves  are  from  Mie  theory, 

which  is  the  exact  result  for  an  uncoated  sphere  while  the  dashed  curves  are  from 
the  Marston-Kingsbury  Physical  Optics  Approximation  (JOSA,  1981).  Figure  4 
also  shows  results  for  a  thinly  coated  bubble  (the  solid  curve  in  this  case)  where  the 
Mie  theory  (for  an  uncoated  bubble)  is  the  curve  with  short  dashes.  The  fine 
structure  is  larger  in  amplitude  for  the  coated  sphere  than  for  the  uncoated  sphere 
but  in  both  cases  the  coarse  structure  is  well  approximated  by  the  POA  Figure  6 


FIGURE  A2 
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shows  the  j  =  1  (E  perpendicular  to  scattering  plane)  case  for  ka  =  500  and  is  to  be 
compared  with  Fig.  3  for  the  corresponding  j  =  2  case.  It  is  seen  (as  noted  in  my 
various  previous  publications)  that:  (i)  the  fine  structure  is  much  smaller  in 
magnitude  for  j  =  2  than  for  j  =  1;  (ii)  the  Mie  results  averaged  near  0C  =  82.8°  are 
closer  to  the  universal  POA  value  of  lj(0c)  =  1/4  for  the  case  j  =  2  than  for  j  =  1;  (iii) 
for  both  j  =  1  and  j  =  2  the  underlying  coarse  structure  is  well  approximated  by  the 
POA,  and  (iv)  both  the  Mie  theory  and  the  POA  show  only  small  differences  in  the 
locations  of  the  coarse  peaks  between  the  cases  j  =  1  and  j  =  2  so  that  it  is  unlikely 
that  these  difference  a  would  be  useful  measuring  ka. 

For  light  from  a  He-Ne  laser  X^T  =  633  nm,  =  (3/4)  *air>  811(1  k  = 
2x/Xwater  Th®  relation  between  the  ka  and  the  bubble  radius  a  and  diameter  2a  is 
listed  below  in  Table  I. 


Table  I 


ka 

radius  =  a  (pm) 

diameter  =  2a  (pm) 

100 

7 JS 

15 

500 

38 

76 

1000 

75 

150 

2500 

189 

378 

3.  Difficulties  Measuring  the  Size  of  Objects  from  Limited  and  Relative  Scattering 

Data  at  Large  ha- 

Inspection  of  Figs.  2-6  show  that  near  0C  the  coarse  structure  is  one  or  more  order(s) 
of  magnitude  broader  than  the  fine  structure.  The  fine  structure  spacing  A0gn#  at 
82.8°  is  approximately  (Langley,  Marston,  Applied  Optics  1984)  for  this 
wavelength: 


A0gM  (degrees)  ■  (1/45)  (a  mm)'*, 


(2) 


Also  our  Mie  calculation  near  forward  scattering  show  the  spacing  of  the  forward 
diffraction  peak  is  of  the  same  magnitude  as  A8gne  given  by  this  approximation. 


Hence  because  of  the  large  size  of  the  bubbles  to  be  studied  and  because  very  closely 
spaced  detectors  are  not  feasible  with  the  Wvatt  instrument,  there  is  no  point  in 
trying  to  use  either  the  fine  structure  or  the  spacing  of  the  forward  diffraction  peaks 
to  measure  the  sizes.  At  best,  we  can  hope  to  minimize  the  complications  due  to  fine 
structure  by  taking  j  =  2  and,  for  each  detector,  by  averaging  over  a  small  range  of 
scattering  angles  (see  item  6  of  letter  of  April  25).  The  point  is  that  even  if  the 
detectors  could  easily  resolve  the  fine  structure,  for  the  purpose  inverting  scatt 
data  it  would  be  impossible  to  sav  which  fine  fringe  was  being  viewed!  This  is 

because  physical  considerations  limit  the  spacing  between  detectors  to  about  5° 
which  is  »  A0fin#. 


Inversion  Procedure  Based  on  Coarse  Structure- 


For  the  reasons  noted  above  it  is  fortunate  that  the  coarse  structure  exists  and  that  it 
has  a  well  defined  starting  angle,  0C  =  82.8  deg.  As  outlined  below  it  should  be 
possible  with  a  few  detectors  placed  close  to  0e  to  measure  a  for  the  size  range  of 

interest  roughly  ka  =  100  -  2500,  and  if  necessary  down  to  ka  ■  25. 

Figure  7  shows  the  normalized  irradiance  I2  (0  =  70  deg)  plotted  as  a 


function  of  ka.  The  plot  is  given  by  the  POA  and  hence  does  not  include  any  fine 
structure.  Since  a  perfectly  reflecting  sphere  has  an  irradiance  proportional  to  a2 
the  actual  irradiance  and  detector  output  is  proportional  to  a2  times  the  function 
plotted  in  Fig.  7.  According  to  the  physical  optics  approximation  Ij  (82.8  deg)  *  0.25 

so  that  the  ratio  of  the  detector  outputs  becomes: 


Wj(0)  ■  Vj(0)/V2(82.8  deg)  -  4  lj(0). 


(3) 


Suppose  v2  (70  deg)  and  v2  (82.8  deg)  are  measured  for  a  given  event  such  that  W2 
(70  deg)  =  4.0  so  that  Eq.  (3)  gives  I2  (70  deg)  =  1.0.  Then,  except  for  the  unusual 

^ se  of  bubbles  outside  the  range  ka  =  100  -  2000,  we  would  expect  the  ka  for  this 
bubble  to  be  given  by  one  of  the  roots  labeled  Aj,  A2,  A3,  and  A4  of  the  condition  I2 

(70  deg)  =  1.0.  Ideally,  to  determine  which  root  we  could  use  the  absolute  magnitude 
of  v2  (82.8  deg);  however,  since  we  do  not  know  the  bubbles  impact  parameter 
relative  to  the  center  of  the  Gaussian  beam  v2  (82.8  deg)  is  of  limited  use  (see  Sec.  B 
of  this  report).  To  discriminate  which  A„  is  appropriate  for  this  event  we  may  use 
the  vj(0)  for  other  detectors  and  apply  Eq.  (3). 

With  this  task  in  mind,  Figs.  8-10  show  I2(0)  for  9  =  75,  65,  60,  and  50  deg  * 
according  to  the  POA.  It  is  evident  that  v2  (75  deg),  and  hence  I2  (75  deg),  would  be 
especially  useful  in  this  task  except  when  the  experimental  I2  (75  deg)  lies  close  to 
one  of  extremum  of  Fig.  8.  Thus  the  roots  A3  and  A4  may  be  clearly  distinguished 
from  Aj  and  A2  by  using  the  measured  value  of  I2.  By  using  a  measured  value, 
via  Eq.  (3),  of  either  I2  (65  deg)  or  I2  (60  deg)  one  can  say  with  certainty  which  of 
A}  or  A2  is  appropriate. 

One  complication  of  this  procedure  is  for  the  case  where  the  measured 
I2  (70  deg)  lies  close  to  one  of  the  extremum  of  Fig.  7.  This  is  the  case  shown  with  roots 
B}  -  B5.  Because  of  uncertainties  in  both  the  measurement  and  the  theory,  a  wide 

range  of  ka  roots  are  now  possible  as  indicated  by  the  dark  band  on  the  ka  axis  below 
roots  B3  and  B4  in  Fig.  7.  Figure  8  suggests  that  I2  (75  deg)  may  be  useful  for 

distinguishing  between  various  roots.  Inspection  of  Figs.  10  and  11  show  that  neither 
measurement  of  I2  (60  deg)  nor  measurement  of  I2  (50  deg)  would  be  of  any  use  in 
discriminating  between  roots  B3  and  B4.  This  is  because  the  band  of  pos-  sible 
values  for  ka  spans  one  or  more  cycles  of  the  coarse  structure  present  in  I2  (when 
plotted  as  a  function  of  ka).  The  ka  spacing  between  these  oscillations  decreases  as  0 

is  reduced.  There  is  a  corresponding  decrease  in  the  angular  spacing  at  fixed  ka  in 

★Fig.ll  for  50  deg  is  omitted  from  this  APPENDIX.  It  is  like  Fig.  10  except 
the  coarse  structure  is  even  more  closely  spaed  in  ka. 


THETA  =  60 
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the  coarse  structure  as  6  decreases.  See  Figs.  3-6.  It  is  dif-  ficult  to  use  measure¬ 
ments  of  vj(0)  (and  hence  lj(0))  for  0  £  60  deg  since  it  will  be  difficult  to  say  with 

certainty  which  of  the  coarse  oscillations  is  being  viewed.  (Recall  that  the  situation 
for  using  the  fine  structure  with  discretely  spaced  angular  detectors  is  even  more 
hopeless.)  Hence  to  make  optimum  use  of  the  coarse  struc-  ture  I  recommend  that  as 
many  detectors  as  possible  be  spaced  between  9  =  65  deg  and  9  =  82.8  deg.  To  reduce 
the  complications  of  averaging  over  the  fine  fine  structure,  these  should  be  placed  on 
the  great  circle  which  is  contained  in  the  plane  of  the  E  polarization  of  the  incident 
wave.  This  corresponds  to  i  =  2  scattering.  This  inversion  procedure  mav  be  aug¬ 
mented  with  proper  use  of  the  absolute  amplitude.  See  Section  B2  of  this  report 

Summary  of  Recommendations  on  Detector  Placement  - 

In  addition  to  the  detectors  placed  as  noted  in  item  4(aMd)  of  my  letter  of  April  25, 1 
recommend  that  as  many  as  possible  of  the  following  detectors  be  placed  in  the  j  =  2 
great  circle.  These  are  listed  in  roughly  the  order  of  priority. 

Table  II.  Recommended  Detectors 


I.D.  No. 

Scattering  angle  at 
the  center  of  the 
detector  (deg) 

Full  angular  width 
viewed  by  the 
detector  (deg) 

1. 

82.8 

1.8  to  3 

2. 

(-)82.8* 

1.8  to  3 

3. 

75 

1.8 

4. 

70 

1.8 

5. 

78 

1.8  to  3 

6. 

72.5 

1.8 

7. 

65 

1.8 

8. 

65 

-1 

"Detector  2  placed  in  opposite  half  arc  as  detector  1.  The  other  detectors 
may  be  in  either  half-arc  in  the  plane  which  contains  the  E  field  of  the 
incident  beam. 
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t: 


r.j 

TV 


Note  that  the  angular  width  cannot  be  large  for  the  detectors  placed  with  0  $  75  deg 
since  otherwise  they  will  not  resolve  the  coarse  features  for  the  larger  of  the 
anticipated  bubble  size.  It  is  assumed  that  these  detectors  will  not  all  be  needed  for 
the  smaller  of  the  bubbles  sizes  so  that  it  is  not  important  to  average  over  the  fine 
structure  of  a  small  bubble  for  detectors  with  0  £  70  deg.  All  of  detectors  5-8  are  not 
essential  if  detectors  1-4  are  available. 


r3 


KH 


K3 


6.  Corrections  to  the  PQA- 

The  POA  prediction  that  the  average  value  I2  (82.8  deg) »  0.25  begins  to  break  down 
for  very  large  ka.  This  is  evident  by  inspection  of  Fig.  5(a).  Hence  it  will  be 
necessary  to  include  some  slowly  varying  function  of  ka  as  a  correction  factor  in 
the  right-hand-side  of  Eq.  (3).  This  factor  may  be  determined  from  Mie  theory. 

The  physical  reasons  for  this  breakdown  at  large  ka  are  discussed  in  Langley- 
Marston  (Applied  Optics,  1984).  The  aforementioned  correction  factor  should  not 
significantly  affect  ability  to  use  either  the  procedure  outlined  above  nor  the 
modified  procedure  listed  below. 

Furthermore,  because  of  the  finite  angular  width  of  the  detectors  and 
approximations  intrinsic  to  the  POA.  the  oscillations  of  the  corresponding  detector 
outputs  will  be  somewhat  less  than  those  displayed  in  Fig.  7-11.  This  reduction  will 
be  moat  significant  at  the  larger  values  of  ka  and  the  smaller  values  of  0  since  the 
half-width  of  the  coarse  angular  quasiperiod  becomes  similar  in  magnitude  to  the 
detector  angular  width.  This  is  evident  from  inspection  of  Fig.  3-5.  In  a  first 
approximation,  it  should  be  possible  to  account  for  this  reduction  by  integrating  the 
Mie  result  over  the  full  angular  width  of  the  aperature;  however,  this  procedure 
assumes  the  effects  of  vignetting  may  be  neglected.  See  item  6  of  my  letter  of  April 
25.  For  a  discussion  of  the  significance  of  vignetting  on  the  measurement  of  an 
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y. 


k. 


object's  angular  spectrum  (in  this  case  the  scattering  pattern  of  a  bubble)  see  J.  W. 
Goodman,  Introduction  to  Fourier  Optics  (McGraw-Hill,  1968)  Figure  5-6  et.  seq. 

B.  The  Value  and  Use  of  the  Absolute  Magnitude  of  the  Voltage  Output  of  the  Detector 
Placed  at  82.8° 

1.  Introduction- 

In  general,  the  voltage  output  of  a  detector  should  be  of  the  form 

Vj(0,t)  =  C(0j)  expt-r^tyb^IjWa2  (4) 


P; 


where  b  is  a  beam  radius  parameter  for  the  Gaussian  beam,  r(t)  is  the  time 
dependent  distance  of  the  center  of  the  bubble  from  the  axis  of  the  Gaussian  beam,  Ij 

is  the  normalized  irradiance  as  discussed  previously,  and  C(0j)  (which  is 
proportional  to  the  laser  power)  is  determined  during  system  calibration.  It  is 
assumed  that  the  bubble  radius  a  is  much  smaller  than  the  hcam  radius  b.  Fora 
given  event  the  vj  takes  on  their  maximum  value  when  r(t)  is  minimized;  call 

this  minimum  r(t)  the  impact  parameter  for  this  event  y.  Hence  y  is  the  smallest 
value  of  the  distance  of  the  center  of  the  bubble  from  the  center  of  the  beam  as  the 
bubble  drifts  through  the  beam.  Note  that  as  in  Fig.  1  rit^^)  =  y.  The  problem  is 

that  y  is  not  known! 

As  shown  below  the  value  of  vj(0,  t^^)  ■  vj(0)  may  be  used  to  simplify  the 


inversion  procedure  proposed  in  Sec.  A.  4  of  this  Report.  It  may  also  be  used  to  check 
the  final  size  distribution. 


2.  Simplified  Inversion  Procedure  Based  on  the  Value  of  vo  (82.8  deg)~ 
According  to  the  POA,  I2  (82.8  deg)  =  0.25  so  that  Eq.  (4)  gives 


[V2<82.8  deg)]/0.25  C  «  (a^)2  =  exp(-y2/b2)  a2, 


(5) 


Note  that  if  y  =  0,  a  =  &min  but  that  generally  a  2  a^  so  that  from  the  experimental 
value  for  v2  (and  calibration  data)  a  minimum  size  for  the  bubble  mav  be 
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determined.  This  value  should  be  computed  for  each  event  used  to  eliminate  some 
of  the  possible  roots  illustrated  in  Fig.  7.  Thus,  for  example,  we  could  rule  out  the 
roots  labeled  Aj  and  A2  if  ka^,,  £  900  so  that  we  need  only  determine  which  of  A3 
or  A4  is  the  appropriate  rot  from  data  on  v2(0)  for  other  9.  Similar  arguments 
could  be  applied  to  case  B  shown  in  Fig.  7.  These  decision  making  processes  could 
be  done  by  a  computer  based  algorithm. 

3.  Use  of  v2  (82.8  deg)  and  the  Empirical  Bubble  Radius  a  to  Define  a  Sample 
Volume- 

The  procedure  outlined  in  Sec.  A4  and  B2  may  be  applied  on  an  event-by-event 
basis  to  determine  the  relative  concentrations  of  bubbles  of  different  radii  within 
some  radius  interval,  for  example,  Aa  =  10  pm.  Call  this  concentration  n(a)  = 
bubbles  per  Aa.  Let  T  be  the  time  duration  for  the  accumulation  of  this  set  of  event 
data,  for  example,  we  may  take  T  =  60  sec.  Suppose,  however,  it  is  desired  to  know 
the  absolute  bubble  concentration  N(a)  =  bubbles  per  Aa  per  volume  of  water.  It  is 
evident  that 

N(a)  =  n(aW,  (6) 

where  V  is  the  volume  of  water  swept  through  the  viewed  area  during  the  time 
interval  T  of  this  data  sample.  During  this  time  interval  we  suppose  that  the 
average  water  velocity  u,  measured  relative  to  the  detector  housing,  is  known.  Let 
ymBX  denote  the  maximum  impact  parameter  for  a  detected  bubble.  This  number 

influences  the  viewed  volume  as  noted  below.  Let  L  denote  the  length  of  beam 
viewed  by  the  detector  located  at  82.8  deg.  It  will  be  assumed  that  L  »  ymax  and 

L  »  b  where  b  is  the  Gaussian  beam  parameter  as  in  Eq.  (5).  Then  the  volume  of 
water  viewed  during  this  time  interval  is  approximately 


V  -  (2  ymax  L)  uT  I  sinyi, 


(7) 
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where  y  is  the  angle  between  u  and  the  axis  of  the  laser  beam.  For  example,  if  the 
background  flow  velocity  is  perpendicular  to  the  laser  beam  y  =  90  deg,  Eq.  (7)  may 
also  contain  small  geometric  correction  factors  which  depend  on  (ymaj/L)  and  on 

sin(82.8  deg). 

One  difficulty  with  the  spectrometer  of  the  Wyatt  Technology  design  (for  the 
present  application)  is  that  the  viewed  area  (2  ymaz  L)  may  not  be  well  defined. 

This  difficulty  was  avoided  by  Gowing  and  Ling  by  use  of  an  optical  mask  to  define 
the  sample  area.  Others,  including  myself,  have  proposed  the  use  of  coincidence 
techniques.  Unfortunately,  none  of  these  techniques  are  easy  to  put  into  effect 
I  will  assume  that  L  is  known  (via  Wyatt's  optics)  and  that  the  problem 
reduces  to  determining  the  effective  beam  half-width  ymu.  Inspection  of  Eq.  (5) 

shows  that  the  impact  parameter  for  a  given  event  is 

y  =  W2  ln(a/amin)]1/2.  (8) 

We  mav  set  ymax  during  the  data  analysis  bv  discarding  all  events  for  which  y  > 
ymax.  I  assume  that  the  electronics  is  sufficiently  sensitive,  and  the  background 

noise  sufficiently  low,  that  the  smallest  bubble  of  interest  will  trigger  an  event 
provided  y  <1  ymax.  As  a  check  on  the  data  analysis  procedure,  the  y  values  for  the 

accepted  events  should  be  uniformly  (and  randomly)  distributed  in  the  interval  0  <, 

y  ^  ymax- 

4.  The  Importance  of  the  Laser  Power- 

I  must  emphasize  again  that  if  the  procedures  outlined  in  Sec.  B2  and  B3  are  to  be 
used  it  is  essential  that  (i)  proper  calibration  data  be  obtained,  and  (ii)  any 
fluctuations  in  the  beam  power  be  taken  into  account  by  uniformly  adjusting  the 
C(0  j).  I  recommend  that  either  laser  output  power  be  monitored  or  the  power  in  the 
beam  dump  be  measured. 


m-Tm 


I  have  also  outlined  an  operational  check  based  on  the  distribution  of  the  measured 
v2  (82.8  deg)  as  well  as  the  inferred  distribution  of  the  radii  a  for  the  sample  of 

events  under  consideration.  This  procedure  assumes  that  the  impact  parameters  y 
are  randomly  distributed  in  the  interval  0  £  y  £  ymax.  The  procedure  makes  use  of 
an  effective  response  function  which  corresponds  to  the  distribution  of  a^^  from 

Eq.  (5),  which  you  would  have  if  the  bnhhles  were  all  of  the  same  size.  This 
distribution  function  gfag^jj/a)  may  be  calculat  d  from  first  principles  if  ymax  is 
known.  Let  the  number  density  for  the  distribution  of  the  (in  a  sample  with 
random  values  of  a)  be  given  by  nfa,^).  Then  we  expect  that  ntamJ  should  be 

given  by  a  convolution  of  n(a)  with  g.  Note  that  in  the  special  case  of  a 
monodispersion  of  radius  bq,  n(a) «*  8(a  -  ao)  where  5  is  a  Dirac  8  function.  Then 
the  convolution  gives  nfa,,^,,) «  g(amin/ao)  as  expected.  In  principle  it  mav  be 
Maaiblfl  to  determine  n(a)  from.  by  a  deconvolution  procedure  but  the 
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